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Abstract 

A review is presented of past and recent attempts to build multiquark states within current 
models already describing ordinary mesons and baryons. This includes: coherence in the chro- 
momagnetic interaction, tetraquarks with two heavy quarks, Steiner-tree models of confinement, 
and hadronic molecules, in particular in the hidden-charm sector Some emphasis is put on the 
difficulties encountered when extrapolating toward higher configurations the dynamics of confin- 
ing forces, starting from the simple case of a quark and an antiquark, or three quarks in a colour 
singlet. 

1 Introduction 

In principle, especially in front of this audience of the International Few-Body Conference, the task 
of describing multiquark spectroscopy should be rather easy: a model is first tuned to fit the spectrum 
and properties of mesons and baryons, and then applied to multiquark configurations. The main effort 
would hence be on the subtleties of the few-quark problems, which might differ from other few-body 
systems by the confining character of the interaction. 

However, the procedure is not that straightforward, because the quark dynamics governing mesons 
and baryons is not very well known, and, within any assumed model, the extrapolation towards mul- 
tiquarks involves severe uncertainties. 

This talk is organised as follows. In Sec. [2l some milestones in the multiquark history are re- 
minded. In Sec. [3l binding mechanisms are discussed in the framework of quark dynamics, with 
emphasis on spin-independent forces. In Sec. |4l the Steiner-tree model of confinement is applied to 
multiquarks. Section[5]is devoted to some conclusions. 

Due to space limitations, the list of references will necessarily be limited. Apologies are presented 
to the many colleagues whose interesting work is implicitly taken into account, but not explicitly cited. 
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2 History 



2.1 Brief survey 

The search for exotic hadrons started even before the quark model, though the "exotic character" was 
not very precisely defined. For instance, it was observed that hyperons with baryon number B = +1 
have a negative strangeness number S, once 5 = +1 is arbitrarily attributed to K^. Some indications 
came, however, on a baryon with S = +1, the so-called "Z" resonance. See, e.g., early issues of 
Review of particle properties, cited in the last one [1]. Eventually, the analysis of the KN data did not 
confirm the existence of the Z baryon. Regularly, indications have been reported about "dibaryon" 
resonances, but, again, no firm confirmation was ever reached. 

Today, by "exotic", one means a state whose quantum numbers and main properties cannot be 
matched by a quark-antiquark [qq] or a three-quark (qqq) configuration. In explicit quark models, 
and in the data, the indications are often in a sector of crypto-exotic states, with the same quantum 
numbers as for ordinary hadrons, but intriguing properties. For instance, a four-quark state with 
isospin / = and spin-parity = 0+ competes with an orbital excitation of {qq) where a spin 
triplet couples to an orbital momentum L = 1, to form = 0^. The problem is whether there 
is evidence for supernumerary states. Otherwise, the dominant {qq) structure is simply modified by 
higher Fock components. 

The experimental situation had unfortunately ups and downs: baryonium, dibaryon, pentaquark 
or molecule candidates have been reported, and usually not confirmed. Some scepticism tends thus to 
prevail among non-experts. On the theory side, also, the situation is seen by outsiders as rather con- 
fuse, with promising scenarios for multiquark binding often not confirmed by more detailed studies. 

Diverse mechanisms have, indeed, been proposed for building hadrons beyond {qq) or {qqq): 
chromomagnetism, chiral dynamics, favourable symmetry breaking, etc., but none of them has been 
fully successful and convincing, so far. There is also a recurrent problem attached to almost any idea 
about multiquark binding: once a state is proposed, with due caution, followers jump on the idea and 
list dozens of "states", most of them lying well above their lowest dissociation threshold, and thus 
likely falling apart immediately and not visible as peaks. 

For example, beautiful mass formulas have been written in the late 70's to describe hadrons made 
of clusters whose mass is shifted by colour-dependent hyperfine forces. These formulas were elab- 
orated using sophisticated group-theoretical techniques. However, estimating hyperfine effects does 
not address the problem of resisting spontaneous fall-apart into two colour singlets. Hence any mul- 
tiquark above the dissociation threshold should be considered with a lot of caution. Jaffe and Low 
|[2l made this point more rigorous. Multiquark states with an artificial confinement constitutes a basis 
of "primitives" for describing the short-range hadron-hadron dynamics, and usually have little to do 
with candidates for new states in the spectrum. 

2.2 Exchange of particles 

We refer here to the review by Roy [3]. The conventional approach to the hadron-hadron interaction 
was based on the Yukawa mechanism at low energy, and Regge poles at higher energy, i.e., the ex- 
change of particles in the t-channel. The hadronic reactions without actual particles to be exchanged 
are thus expected to be suppressed. This is the case, for instance, of backward tt+tt" elastic scattering, 
seen as a forward vr+vr^ — )• 7r^7r+ reaction, which involves charge Q = 2 exchange. 

On the other hand, hadronic reactions can be described by their s-channel content. The interaction 
is stronger in channels with resonances than without resonances. The consistency between these 
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Figure 1: s- and t-channel content for some hadronic reactions: (a) meson-meson, (b) meson-baryon, 
(c) and (d) baryon-antibaryon, and (e) baryonium-baryon scattering. 

two pictures was named duality. The reactions with large cross sections must have actual particles 
exchanged in the i-channel and resonances in the s-channel. This concept was initiated before the 
quark model, but became more convincing when duality diagrams were drawn in term of quark lines, 
as illustrated in Fig. [T] 

When considering the case of nucleon-antinucleon scattering, to which ordinary mesons con- 
tribute in the t-channel, Rosner pointed out the need for new type of mesons, preferentially coupled to 
baryon-antibaryon, as their partners in the s-channel, see Fig.[T] Later on, dynamical models were in- 
vented to tentatively describe such "baryonium" states, as seen below. Note that if you now enforces 
an ordinary meson in the s-channel, you need a baryonium in the t-channel to satisfy the duality 
requirements. 

As pointed out in [3 ], duality does not escape the Pandora-box syndrome that is encountered in 
most approaches to multiquarks. Once baryonium is accepted as a new hadron, the pentaquark au- 
tomatically emerges. Consider, indeed, baryonium-baryon scattering, with ordinary mesons (quark- 
antiquark) exchanged in the t-channel. Then the partner s-channel consists of pentaquark, as seen 
in Fig. [T] Even more complicated structures appear if, in turn, a pentaquark enters a scattering with 
another hadron. 

2.3 Baryonium 

In the 70s, the baryonium states anticipated by duality arguments were tentatively seen as peaks in the 
antiproton cross section (M > 2 nip), or peaks in the annihilation at rest pp ^ ^ + X (Mx < 2 nip) 
m. They were not confirmed, unfortunately, in later experiments with higher statistics and better 
resolution, in particular at the LEAR facility of CERN. 

A somewhat simple picture of baryonium is that of a deuteron-like A^A'^ bound state or resonance, 
benefiting from the attractive potential mediated by the exchange of mesons (5). The authors, however, 
never agreed on the influence of the very strong annihilation on the spectrum. 

Another interpretation was named diquonium or quark-model baryonium, with two localised 
quarks and two locahzed antiquarks separated by an orbital barrier. This [{qq) — {qq)] structure 
explains why the state is preferentially coupled to baryon-antibaryon, rather than to mesons ||6l|4l. 
Such a clustering is, however, assumed and not demonstrated from any dynamical calculation. 

A further speculation envisaged the possibility of a colour sextet for the diquark, this leading to an 
even narrower [{qq) — {qq)]- This "colour chemistry" was also extended to the baryon number B = 1 



sector, with [{qq — (qqq)] states Q. 

2.4 Scalar mesons 

Within models available in the 70s (simple quark models, MIT bag), it was discovered that four quarks 
in S'-wave, experiencing a favourable spin-spin interaction, have a mass which is comparable, or even 
lower, than an orbital excitation of (qq) with a spin triplet S = 1 and an orbital momentum L = 1 
coupled to J = lEm. 

At about the same time, KK molecules were suggested ifTOllTTIl . also candidates for scalar mesons 
near 1 GeV. This was the beginning of an intense activity, with mixing between these various (non- 
orthogonal!) configurations, and further complications due to glueball and hybrid states entering the 
game. 

The main question is whether counting the states in the spectrum is compatible with pure (qq) with 
higher Fock state corrections, or does require extra states. For instance, it is argued in [12] that the 
present data in the meson sector do not call for hybrid, glueballs or molecules, except as corrections 
to a main {qq) component. 

The scalar sector of light mesons is regularly revisited. In the framework of QCD sum rules, it 
is found that the states contain a large fraction of glueballs (13]. A realistic description of the scalar 
mesons and meson-meson interaction requires the contribution of several channels 1,14.1 . 

2.5 Pentaquarks 

Interesting reviews have been written on the pentaquark episode |[T5l [161 . The story started with 
a nice speculation based of chiral dynamics lUTl . leading Nakano and his colleagues to attempt a 
measurement out of the stream of fashion at that time, and to get a positive signal 11811 . The real 
surprise was the wave of positive evidences elsewhere, including experiments having already the data 
on tape for several years, that nobody had to curiosity to use to search for exotic baryons. However, 
the pentaquark was not confirmed in high-statistics experiments using sophisticated and expensive 
detectors with good particle identification. 

Some studies were also done in lattice QCD, but somewhat hastily, with a variety of conclusion 
about the existence of the pentaquark and about its plausible quantum numbers. Clearly, more time 
would have been required to better distinguish between a genuine bound state or resonance and a 
state simply occurring from the discretisation of the continuum in the finite volume of the lattice. The 
theoretical activity is much slowed down today, but not completely dead, in particular better studies 
are performed on multiquarks within the framework of QCD sum rules and lattice QCD, as seen in 
several contributions to this Conference, where references can be found. 

2.6 Hadronic molecules 

This is a very long history. In the 50s, the A was described as vrA'^ resonance, i.e., a kind of molecule. 
But this picture was not easily extended to other meson and baryon resonances, and when the quark 
model came, it was very welcome as a progress. 

Today, the quark model is considered as too primitive, especially in the light-quark sector, and is 
given up by many authors, who describe most hadrons as molecules, again, but within the improved 
framework of effective theories. Many interesting contributions have been presented at this Few- 
Body Conference. Any new resonance, especially in the hidden-charm sector, is almost immediately 
described as a new meson-meson molecule, and the Pandora box syndrome seemingly strikes again. 
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However, such a blame cannot be addressed to the X(3872), since it was predicted before its 
discovery, on the basis of the pion-exchange dynamics. See, e.g., [19] for references). If the X(3872) 
were not confirmed as a hadron-hadron molecule and eventually interpreted as being mainly a (cc) 
state, we will feel like a hen sitting on a duck egg, and rather surprised when it hatches. We learned 
to be careful with the molecular interpretation of the hidden-charm states. When higher ijj resonances 
were found, theorists were puzzled by the anomalies in the relative decay rates into DD, D*D + cc. 
and D* D* , and a molecular interpretation was suggested [20| (see, also, \\1\D . But the branching 
ratios were later understood from the node structure of the decaying states Il22ll23ll24l . 

As reported at this conference, the mass and the early-measured properties of the X(3872) are 
very indicative of a DD* + cc. structure. However, the branching ratio of X — )■ '>p{2s) + 7 to 
X J/'ip{ls) + 7 is large, a property that is natural for a radially excited P-state of (cc). 

On the theory side, the problem is whether a satisfactory theory of nuclear forces can be built 
outside the nucleon-nucleon case. Is has been pointed out |[25l that due to the D — D* mass differ- 
ence, the pion propagator carries energy transfer together with momentum transfer, and thus that the 
resulting interaction is non local, and perhaps less effective. In fact, the X(3872) in this approach 
looks like a D — D — tt three-body system, with the pion oscillating back and forth between the 
two heavy mesons, and resonating with each of them. Another problem is that there is a hard core 
in the nucleon-nucleon interaction, preventing them to merge their constituents. No such core exists 
for DD* and other meson-meson systems. Hence the short-range dynamics always plays a role. Of 
course, at very low binding energy, the meson-meson wave function is very extended, but the strength 
of the attraction is not given solely by pion exchange, the exchange of heavier mesons (w, . . .) has 
been introduced (see, e.g., lfT9l ). and this can be seen as a parametrisation of the direct interaction 
between the quarks of each flavoured meson. 

2.7 Mixing dynamics 

This is a current trend in the physics of hadrons: if a first model describes satisfactorily some proper- 
ties of a hadron, but fails for others that are accounted for by a second model, a tantalising improve- 
ment of the picture consists of writing the wave function as 

= cosi9 Tpi + sini) ip2 , (1) 

and this can be extended as to include more components. This game has been played endlessly for 
scalar sates involving radial excitations, hidden strangeness, glueballs and hybrids. In the case of the 
X(3872), Tpi could be a DD* + cc molecule, and Tp2 a charmonium cc(2P) with the same quantum 
numbers and a radial excitation. However, the dynamics of coupled channels dictates that if cos t9 i/ji is 
the leading term, the small admixture sin ■ip2 is not very much governed by the diagonal interaction 
in the second channel. Instead it is given by a folding of the leading component and the transition 
operator. In particular, if the latter is nodeless and the former smooth, a node hardly appears in ip2. 

A toy model illustrating this property is an S-wave analogue of the X(3872) system, with a (DD) 
(the difference between D and D* is ignored) channel weakly bound by a Yukawa potential, and a 
(cc) channel with a standard Coulomb-plus-linear potential, such that before mixing, the threshold 
is at 3.8 GeV, (DD) at 3.77 GeV and ^^(25") at 3.57 GeV. The mixing is crudely mimicked by a 
local interaction with a Gaussian shape; The mass of the (DD) is slightly shifted (more details will 
be given elsewhere), and the wave function acquires a (cc) component, displayed in Fig. ??, which is 
nodeless. 

The current belief is that the eigenstates mix unchanged, dominantly by affinity of neighbouring 
unperturbed energies, once a gate is open between the two channels. This is inspired by the denomi- 



5 




Figure 2: Mixing of {DD) and (cc) in a simple local model. The solid lines represent the reduced 
radial wave function for the main {DD) and the small (cc) components, as obtained from an actual 
coupled-channel calculation. The dashed lines correspond to the naive mixing scheme of Eq. with 
neighbouring unperturbed states, and 9 adjusted to reproduce the same normalisation for (cc). Note 
that the coupled-channel calculation does not produce any node in the (cc) sector. Units are GeV^/^ 
for u{r) and GeV^^ for r. 

nator in the first order correction to a wave function within perturbation theory, namely, in a obvious 
notation, 

(P^-^^ +2^ (0) (0) + 

However, the node structure of the sometimes gives more drastic constraints. Another example 
is 5 — D mixing in ordinary charmonium dynamics. It is almost ever considered that il>{lD) pref- 
erentially mixes with tp{2S) which is very close. But when one actually computes the small 5-wave 
admixture in tp{lD) using an explicit tensor force with suitable regularisation, one finds a nodeless 
radial function for the admixed S-component in the ^(3770). Note that S-D mixing can be calculated 
analytically for the muonium (/i+e") and the approximations can be tested there. In the above model, 
the ground-state has two nodeless components, with the (cc) one dominating. The second state has 
two nodes, one in the main (cc) channel, and another in the small (DD) admixture, as seen in Fig. [3] 
The third one has no nodes. Many variants are possible, if the parameters are modified. 

Back to X(3872), the iIj{2S) to ipi^S) ratio of radiative decays has to be studied within a consis- 
tent picture involving probably several Fock-space components, cc, ccqq, etc. This requires a mod- 
elling of the transition operator creating of annihilating a pair of light quarks, and a few calculations. 

A consistently-managed mixing scheme can give interesting results. An example is the celebrated 
Cornell picture of (cc), supplemented by coupling to real of virtual (including Dg) decay 

channels 112411. It was possible, using this model, to predict before the discovery of the r]c{2s) a 
substantial reduction of the 'ip{2s) — rjc{2s) splitting ll26l . as compared to simple potential models. 
The same effect probably explains why the T — rj^ splitting is slightly larger than estimated by most 
authors: the reduction due the meson-meson channels is weaker here than for J/ip — rjc- Other 
interesting treatments of the configuration mixing can be found, e.g., in [i27ll28ll29l . 
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Figure 3: Reduced radial wave functions for the (cc) (blue) and {DD) (red) components in the above 
coupled-channel model, for the ground-state (dotted lines), first (dashed hnes) and second (solid lines) 
excitations. 

2.8 Diquark clustering 

The concept of diquark is very useful in several branches of particle physics, for instance to analyse the 
baryon-to-meson ratio in multiparticle production. It has been introduced rather early in spectroscopy, 
see, e.g., Ref. |30| for a review. However, some of the pioneering works are sometimes ignored in 
recent rediscoveries of the diquarks. A first problem to which diquarks provide a sufficient solution, 
is why Regge trajectories (squared mass vs. spin J) are linear with the same slope for baryons as 
for mesons. For two-body systems, the linear character is reproduced in many pictures, e.g., H2 = 
^(m^+pf )^/^ + Ari2. Thus the equality of slopes comes automatically if the quark-diquark baryons, 
[q — {qq)], are the partners of the quark-antiquark mesons, [q — q\- 

However, it is not necessary to introduce diquarks by hand to obtain the equality of slopes. In the 
symmetric quark model, the baryon analogue of H2 reads = ^(m^ + pfY^"^ + V3, where V3 is 
A ^ rij/2 or the y-shape potential (discussed in Sec.|4l). The Hamiltonian H^, for the large angular 
momentum, gives linear Regge trajectories for baryons, with the same slope as mesons. There is a 
dynamical clustering, or, say, a spontaneous breaking of symmetry when J increases IIBTI . 



Figure 4: Schematic picture of a rotating baryon at low angular momentum J and higher J. 

The best argument in favour of the quark-diquark model perhaps lies in the problem of missing 
resonances. In a recent review on baryons 1,32 J . very few states are tentatively interpreted with both 




large J 
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degrees of freedom (i.e., Jacobi variables x and y) excited. On the other hand, many states predicted in 
the three-quark model are not observed, and this was the subject of contributions to this Conference. 
The most striking state with double excitation in the three-quark model is 



^! = X xy exp[-a (x^ + y^)/2] , 



(3) 



(here in the case of harmonic interaction, but a state with the same symmetry exists for other confining 
potentials). For experts, it is named the [20, 1+] multiplet. Its fully antisymmetric orbital function, 
with N = 2 degrees of excitation, is associated with an antisymmetric spin-isospin wave function, 
and an antisymmetric colour wave function. It is absent in the quark-diquark picture. Experimentally, 
it is not (yet?) seen. 

If the diquark model is taken seriously, it can be extrapolated outside the framework of baryon 
spectroscopy. We have seen that the late baryonium states have been described as a diquark and 
an antidiquark. More recently, cq or cs diquarks and the associated antidiquarks were introduced to 
describe the X,Y or Z states seen in the hidden-charm spectrum. See, e.g., |[33l . The light pentaquark 
was also described as [{qq) — (qq) — q] l[34l . 

Obviously, the Pandora box syndrome becomes threatening here too. In particular, three diquarks 
can well build a dibaryon. It should be checked whether in models describing the Y mesons as 
(cs) — (cs), the configuration (cs)^ is not below the threshold for (sss) + (ccc), since it would be 
embarrassing to pay the price of a very exotic dibaryon to explain the Y meson. Years ago, a light 
dibaryon (or "demon deuteron" ) was shown to be a consequence of the diquark model |35 ]. HI 

3 Binding mechanisms in simple quark models 

3.1 Chromomagnetic binding 

In the 70s, the chromomagnetic interaction. 



inspired by the Breit-Fermi interaction in atoms, was shown to account for the observed hyperfine 
splittings in ordinary hadrons, including the S — A mass difference. There is an analogue in the 
bag model. In fBll, Jaffe pointed out the H(uuddss) in its ground state has a chromomagnetic 
interaction whose cumulated (attractive) strength is larger than for the lowest threshold K{uds) + 
A(uds). Note the non-trivial colour-spin algebra. In a model with pure spin-spin, Vss oc X] ^i-^j^ ^ 
J = tetraquark receives a cumulated strength = —6, whereas each scalar meson of its 

threshold experiences Ylii<j = ~3, and thus there is no obvious excess of attraction. At best a 
small gain from a favourable spatial distortion. Remember also the positronium molecule, governed 
by V = J2i<j dij/'^ij'- its binding is not obvious because the cumulated strength is 9ij = —2 for 
the molecule, and ^ gij = —2 also for the threshold made of two atoms. The miracle in (01 is that 



If the quark masses rrii are identical (SU(3)f symmetry) and if the short-range correlation factors 
{6^^^ {vij)) are assumed to be the same in multiquarks as in ordinary baryons, then the H is predicted 

'There is an error in 1351 about the quantum numbers of the orbital wave function, but this does not remove the issue of 
dibaryons in diquark models. 




(4) 




(5) 
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about 150 MeV below the threshold. The H was searched for in more than 15 experiments, with 
negative results, including double-A hypernuclei. 

In 1987, Lipkin, and independently Gignoux et al. 1371 l38l |391 pointed out that the same mecha- 
nism would bind an exotic pentaquark (the word was invented in this context) 

PiQqqqq) < (Qq) + {qqq) , (6) 

where Q = c, b and (qqqq) is a combination of u, d, s, in a SU(3)f flavour triplet. The same binding of 
— 150 MeV is found if Q becomes infinitely heavy. A search for the P r in an experiment at Fermilab 
turned out not conclusive [40, 41]. Silvestre-Brac and Leandri listed other candidates B2ll . 

However, further studies of the chromomagnetic model fl3l l44l l45l |46]| indicated that the correc- 
tions tend to moderate and even to cancel the binding effect due to the group-theoretical result ([5]). In 
particular: 

• There is more kinetic energy in a dibaryon than in two baryons (5 Jacobi variables vs. 4), 

• SU(3)f breaking is not favourable. The A resists much better than H. 

• The matrix elements {5^^'^ {uj)) appreciably smaller in multiquarks than in ordinary hadrons. 

Nevertheless, the simple chromomagnetic Hamiltonian (HJl might still reveal some surprises. For 
instance, if diagonalised with realistic values for the flavour dependent strength factors {6{fij)/ {niimj), 
it provides the {ccqq) lowest configuration with J^^ = 1++ the main features required to describe 
the X(3872) Ii47ii48il . 

3.2 Chromoelectric binding 

We adopt here the language of potential models but we believe that the results are much more general. 
The main advantage of potential models is the possibility of switching on or off some contributions to 
single out the most effective one for binding. 

A remarkable property of the spin-independent interaction among quarks i?, flavour independence, 
which induces interesting symmetry breaking effects. 

Remember that symmetry breaking tends to lower the ground-state energy. The simplest example 
is /i = + + Xx in one-dimensional quantum mechanics, with a ground state energy e(A) = 
1 — which is always below e(0). But this is very general. If 

H{\) = i^even + A Hodd , (7) 

the variational principle applied to -ff(O) with the even ground state ^'(0) of -ffeven as trial wave 
function gives 

E{\) < (^(0)|i/(A)|^(0)) = S(0) . (8) 

One can apply this result to few-body systems with a variety of symmetries for which the labels "odd" 
and "even" make sense, in particular particle identity and charge conjugation. But stability is a com- 
petition between a configuration with collective binding and another configuration where the system 
is split into separate decay products. The threshold also benefits from the symmetry breaking, and 
very often more! In this case, stability deteriorates, though the energy of the compound configuration 
decreases. 

For instance, consider the barely bound (e"*", e"*", e~, e~) molecule, or any rescaled version with 
the electron mass replaced by another mass ji, and move to configurations involving two different 
masses. Then it is observed, and proved, that: 
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• for (Af+, M~, m~): binding deteriorates and is lost for M/m > 2.2 B9l . 

• for {AI'^ , , m~ , m~): the binding is improved. 

But the Coulomb character matters little. What is important, is that the potential does not change when 
the masses are varied. Hence, a similar behaviour is observed for any four-body problem with flavour 
independence. The splitting (the potential V is assumed to be symmetric under charge conjugation, 
and independent of the masses) 



1 1 

4M 4m 



(Pi - P2 + P3 - Pa) 



implies for the ground state 

E{M,M,fh,fh) < E{fi,n,fi,fi) , 2^^^ = M"^ + m"^ , (10) 

but, meanwhile, the threshold energy remains constant, as (M, fh) and (/i, /i) have the same reduced 
mass. So the stability is improved. 

Explicit quark-model calculations have been carried out to illustrate how this favourable symmetry 
breaking works with flavour-independent potentials. The corresponding four-body problem is rather 
delicate, as most other four-body problems. Remember that when Wheeler proposed in 1945 (the 
paper paper was published somewhat later ll50l ) that the positronium molecule might be stable, a first 
numerical investigation by Ore lISTl concluded that the system is likely unstable, but the following 
year, Hylleraas and the same Ore published a beautiful analytic proof of the stability f52l. 

In current quark models, the main conclusion is that binding a doubly-flavoured tetraquark re- 
quires a large mass ratio, usually (bbqq) or {bcqq). However, a more sophisticated calculation by Jane 
and Rosina ll53l found (ccqq) barely bound. See, e.g., |[54l for a detailed survey of the situation. 



4 Steiner-tree model of confinement 

It should be acknowledged, however, that these early constituent-model calculations suffer from a 
basic ambiguity: how to extrapolate from mesons towards multiquarks. The usual recipe is 

y = -^E^-^-(^^^)^ (11) 

i<j 

which is presumably justified for the short-range part, but not for the long-range part, except for very 
peculiar spatial configurations. In (fTTI) . the normalisation is such that v{r) is the quark-antiquark in- 
teraction in ordinary mesons. Strictly speaking, ([TT]) holds for a pairwise interaction with colour-octet 
exchange. Clearly colour-singlet exchange cannot contribute to confinement, otherwise everything 
would be confined together, but colour-singlet exchange can contribute to short-range terms. More- 
over, there aie very likely three-body and multi-body forces in baryons and multiquarks. 

In the case of baryons, it was suggested very early ||55l|56l that the potential generalising the Unear 
confinement of mesons is 

3 

Vy = CF min ^ r^j . (12) 

1=1 
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Figure 5: Confinement of mesons and baryons. Tiie minimum over the quark permutations gives the 
flip-flop potential. 

This was often rediscovered in the context of models (adiabatic bag, flux tubes), or in studies dealing 
with the strong-coupling regime of QCD [57|. Estimating the baryon energies and properties with 
the potential (fT2l ) is a very interesting 3-body problem. However, the result differ little from these 
obtained using the colour-additive rule, which for baryons reduces to the "1/2" rule, namely 

a 

^3 = ^ (n2 + ?'23 + ?'3l) • (13) 

This y-shape interaction has been generalised to tetraquarks. At first, this looks as an astute guess, 
but it was later endorsed by detailed lattice QCD ||58l . including the interplay between flip-flop and 
connected Steiner tree. See, also, ||59l for a study within AdS/QCD. The confining potential reads 

V = min{(fi3 + ^24,^14 + c?23, V4} , (14) 

14 = mill (lluiSill + ||u2Si|| + |jsiS2|| + |[S2^^3|| + ||'S2t'4|| ) , 

corresponding schematically to the flux tubes in Fig. [6] 

V\ V3 Vl 



#^ ^ 

V2 V4 V2 



Figure 6: Generalisation of the linear potential of mesons to tetraquarks: the minimum is taken of the 
flip-flop (left) and Steiner tree (right) configurations. 

A first study of the tetraquark spectrum with this potential concluded to the "absence of exotics" 
II6OI . but a re-analysis by Vijande et al. [61.1 with a better wave function, indicated that this potential, if 
alone, and free of constraints due to the Pauli principle, gives stability for the equal-mass case (qqqq), 
and improved stability for the flavour asymmetric (QQqq)- It remains to analyse how this stability 
survives antisymmetrisation, short-range terms in the potential, relativistic effects and spin-dependent 
corrections. This is however, very encouraging. 

Two developments came. First, a better understanding (at least by physicists) of the Steiner trees, 
with efficient algorithms |[62l to compute them, and novel inequalities |[63l . For instance, in a Steiner 
tree, a set of point can be replaced by its associated torroidal domain. For instance, in the baryon 
y-problem, the continuous minimisation over the location of the Steiner point s can be replaced by a 
discrete maximisation: the desired Y-shape length, \\svi\\ + \\sv2\\ + Wsv^W in Fig. |7]is the largest of 
the two distances 1 1 sws 1 1 and | i sts 1 1 , where the points and 1 3 , named the Melznak points or torroidal 
domain of {vi,V2}, form an equilateral triangle with vi and V2- 
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Figure 7: The continuous minimisation over the location of s is replaced by the discrete maximum 
max (llswall, ||si3||}, where ^3 and is make an equilateral triangle with the quarks vi and V2- Here, 

V3 = \\sVi\\ + \\SV2\\ + WSV^W = \\V3W3\\. 

This property, and the interesting symmetry restoration (even for an asymmetric triangle, the flux 
tubes from the junction to the quarks make 120° angles in Fig. [5ll are related to the theorem by 
Napoleon (see Fig. [8]): if one starts from an asymmetric triangle viV2Vs and builds an external equilat- 
eral triangle such as viV2Ws along each side, the centres of these auxiliary triangles form an equilateral 
triangle. 

The analogue for a planar tetraquark is 

Vs = \\siVi\\ + IIS1V2II + ||s2t'3|l + ||s2'y4|| = |Ilfl2W^34|| , (15) 

and is illustrated in Fig. |9l the length of Steiner tree linking the quarks vi and V2 to the antiquarks V3 
and V4 via the junctions si and S2 is he distance between the two Melznak points W12 and W34,. 

In space, one gets the problem pictured in Fig. [TOl the auxiliary points W12 and w^i are located 
on a circle of axis V1V2 (resp. v^v^), and estimating the minimal Steiner tree can be shown to be 
equivalent to finding the maximal distance between the two circles [63J. This is a standard problem 
in computer-assisted geometry, as applied, e.g., in the cartoon industry. 

Then you can derive inequalities on the Hamiltonian, and recover rigorously in some limiting 
cases the stability established by numerical methods. 

The disappointing observation, however, is that the dynamics is dominated by the simple flip-flop 
term, and the most interesting connected Steiner tree plays a relatively minor role. 

The exercise can be repeated for the pentaquark, using the linear model with minimal cumulated 
length, 

Vp = min(yff, Vst) , 

Vs = mill [rii + Vy(rj,rfc,f^)] , (16) 

i 

Vst = connected Steiner tree , 

i.e., the configurations in Fig.fTT] 

A simple variational calculation gives stability at least for (qqqqq), and {Qqqqq) and {qqqqQ) 
where m{Q) ^ m{q) ||64ll . The other mass configurations remain to be studied. 
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Figure 10: The length of the Steiner tree is the maximal distance between the quark and antiquark 
Melznak circles. 
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Figure 11: Pentaquark confining potential: flip-flop (up) and connected (down). 



This proliferation of stable states in the minimal-length model becomes embarrassing. Very likely, 
the dibaryon will also be found stable. If one believes into this confinement, this means that the role 
of the neglected effects should be investigated with care, in particular: 

• Coulomb forces, perturbative and non perturbative, 

• antisymmetrisation. 

If the constraints of antisymmetrization turn out the main obstacle to multiquark stability, then exotic 
hadrons have to be searched in configurations with quarks of different flavours. 

5 Conclusions 

The problem of multiquark binding is now addressed very seriously with QCD sum rules. Lattice 
QCD and even AdS/QCD. These ambitious but delicate approaches have first confirmed some results 
that were previously obtained empirically, such as the Steiner-tree structure of the linear term of the 
confining interaction. 

Constituent quark models remain, however, in the forefront of investigations, to detect the most 
interesting configurations. With a proper treatment of relativistic effects, and adequate dynamics, 
quark-model calculations can give a very good account of the spectrum and main static properties, 
as illustrated, e.g., by the Bonn group for mesons and baryons [65]. The case of multiquark is of 
course much more difficult, with the mixing of confined channels and hadron-hadron components in 
the wave function. 

On the experimental side, it is hope that the future collider experiments will devote a reasonable 
amount of time to search for exotics with heavy flavour. As shown by B factories, there is a very good 
potential of discoveries in hadron physics within experiments primarily designed for studying other 
aspects of particle physics. 
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